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Abstract. - Chiral Lagrangian and quark- meson coupling models of hyperon matter are used to 
estimate the maximum mass of neutron stars. Our relativistic calculations include, for the first 
time, both Hartree and Fock contributions in a consistent manner. Being related to the underlying 
quark structure of baryons, these models are considered to be good candidates for describing the 
dense core of neutron stars. Taking account of the known experimental constraints at saturation 
density, the equations of state deduced from these relativistic approaches cannot sustain a neutron 
star with a mass larger than 1.6-1.66 Mq. 



Neutron stars are the most compact stellar objects, and 
within them, particles that arc unstable on earth could be 
stabilized at densities above 10 14 g cm -3 . The equation of 
state (EoS) at these densities is still largely unknown due 
to the poor experimental knowledge on hadron dynamics. 
However, the recent measurements of the millisecond pul- 
sars PSR J1614-2230 extend the maximum observed mass 
from (1.67±0.02) M Q [T] to (1.97±0.04) M Q g]. This im- 
poses new constrains on the EoS of dense matter [3J. A 
major question is whether the EoS including hyperons can 
be ruled out by this new observation. Relativistic mean- 
field [H[5] and Brueckner-Hartrcc-Fock 6 8 models have 
predicted that the effect of hyperons in dense matter is 
to soften the EoS and, therefore, to lower the maximum 
mass. A consistent treatment of the Fock term in a rel- 
ativistic approach is however missing. In this Letter we 
analyze the impact of the Fock term in two effective, rel- 
ativistic meson-exchange models related to the quark na- 
ture of hadrons: the chiral Lagrangian model [3] and the 
quark meson coupling (QMC) model |10) . 

These models are based on a relativistic meson-exchange 
theory of the Walecka type, where the mean-field proper- 
ties of nuclear matter are calculated in the Hartree-Fock 
approximation |llj . Considering scalar s, oj, 8, p and 7r 
mesons as well as nuclcons and hyperons, the Lagrangian 
can be written with the usual notations |10j 
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where the sum over the flavor index / means a summation 
over the baryons: N, A, E, S. The isospin operator I and 
the Gamow-Teller operator Gj, are given in Rcf. (It)] . 

In the chiral Lagrangian model [H] , the meson fields are 
related to the properties of the QCD condensate. Scalar s 
and pseudo-scalar 7r fields are associated with the fluctu- 
ations of the chiral quark condensate [12], related to the 
non- linear "Mexican hat" potential, 



(1) 
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deduced from chiral spontaneous symmetry breaking. In 
this model, the chiral partner A\ of the p meson, being 
too heavy, is ignored and the u> meson is chiral singlet. 
The first chiral models failed to reproduce the saturation 
properties of nuclear matter |12j . A better description 
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Table 1: Parameters of chiral and QMC models and saturation properties: density (no), binding energy (Bo), incompressibility 
(Ko), symmetry energy (a s ), symmetry energy slope (Lo) and in-medium Dirac mass for nucleons (Mj^). QMC700 is given in 

Ref. Eg. 
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m s 


9s 
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9 P 


fop 


c 


no 


Bo 


A'o 


a s 
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( r, — 3 "\ 

(rm ) 


(MeV) 


(MeV) 


(MeV) 


(MeV) 


(MeV) 


(MeV) 




MC1-H 


800 


10.0 


5.06 


4.14 


0.0 


1.44 


0.160 


-15.85 


248 


30.0 


79.4 


0.89 


MC1-HF 


880 


10.0 


7.09 


2.65 


4.2 


2.00 


0.160 


-16.01 


271 


30.3 


66.9 


0.81 


MC2-HF 


860 


10.0 


7.31 


2.65 


1.2 


1.90 


0.160 


-16.05 


276 


30.6 


67.9 


0.80 


MC3-HF 


900 


10.0 


6.33 


2.65 


5.6 


1.92 


0.160 


-15.91 


274 


30.2 


63.1 


0.83 


MC4-HF 


880 


10.0 


6.57 


2.65 


5.6 


1.83 


0.160 


-15.92 


274 


30.5 


64.1 


0.82 


QMC700 


700 


11.94 


10.66 


4.167 


0.0 


0.51 


0.159 


-15.74 


338 


29.9 


53.3 


0.74 


QMC-H 


700 


10.52 


8.39 


4.167 


0.0 


0.40 


0.160 


-15.83 


292 


32.7 


89.5 


0.78 


QMC-HF1 


700 


11.90 


10.62 


4.167 


0.0 


0.51 


0.160 


-15.72 


338 


30.4 


54.4 


0.77 


QMC-HF2 


700 


11.30 


9.31 


2.87 


0.0 


0.46 


0.160 


-15.74 


360 


29.8 


72.2 


0.67 


QMC-HF3 


700 


9.78 


6.84 


3.75 


0.0 


1.00 


0.160 


-16.00 


285 


31.4 


69.2 


0.79 



for the saturation is obtained by introducing the scalar 
polarizability of the nucleon k,ms in Mf(s) p"3lll4j . 

M f (s) = M f +g s w f s s+^K NS w f s s 2 , (3) 

which reflects the effects of confinement as originally pro- 
posed in Ref. [IS] . The weight factors w( and w{ will be 
discussed later. 

A further extension of the model consists in replacing 
k ns in Eq. © by k N s(s) = KnsO-+(2s)/(3U)) where /„• 
is the decay constant of the pion [§] . Another extension is 
related to the p and it tensor interactions, where derivative 
couplings lead to a contact term that can be partially or 
totally suppressed by short-range correlations [9][TT]. We 
consider chiral models at the Hartree approximation (such 
as MC1-H in table [1} as well as in the Hartree-Fock ap- 
proximation (such as MCi-HF in table [IJ . The index i in 
MCi-H(F) models refers to different combinaisons of the 
extensions discussed: («=1), constant kns in Eq. ([3]) and 
totally suppressed contact term; (i=2) kns(s) in Eq. ([3]) 
and totally suppressed contact term; (i=3) constant kns 
and partially suppressed contact term; (i=4) kns{s) and 
partially suppressed contact term. More details on these 
extensions can be found in Ref. [5]. 

In the quark-meson coupling model QMC700 the scalar 
potential ^ is quadratic, V(s) = m 2 s s 2 /2, at variance 
with the chiral Lagrangian model the scalar po- 

larizability of the nucleon kns is kept constant as in 
Eq. ([3]); and only s, ixs and p vector mesons are consid- 
ered [k p = 0). The other extensions presented here-before 
are disregarded. In Ref. [10], self-energies are calculated in 
the Hartree approximation while a non-relativistic approx- 
imation of the Fock contribution to the energy is included. 
In this paper, we present modified versions of the QMC 
model in which both Hartree and Fock contributions are 



treated consistently. 

The parameters of chiral Lagrangian and QMC models 
are given in table [TJ In the chiral model, the scalar cou- 
pling constant is fixed from the sigma linear model to be 
g s = Mjsilj-K = 10. Since the dimensionless parameter 
C = J^kjvs/(2Mjv) estimated from lattice simulation is 
around 1.25 [18] . the chiral model fixes C to be consistent 
with this value but a flexibility is retained for reproduc- 
ing saturation properties as explained below. In the QMC 
model, the dimensionless parameter C is smaller and de- 
duced from the nucleon mass in the bag model. The cou- 
pling constant g p is adjusted to fit the symmetry energy 
in the chiral Hartree model MC1-H and in QMC models 
while in chiral Hartree-Fock models MCi-HF g p is fixed 
to satisfy the vector-dominance model giving g p — 2.65 
and the p-tensor coupling constant n p is adjusted to the 
symmetry energy. The axial coupling constant is fixed to 
be gA = 1-25 and the S coupling is fixed to gs = 1. The 
other parameters (m s , g u and C) are chosen to reproduce 
the saturation properties (no, Bq and Kq) reported on ta- 
ble [T] In QMC model, the couplings g s and g^ are fixed 
to reproduce no and Bq; moreover in the parameteriza- 
tion QMC-HF3 the parameter C is modified to decrease 
the value of Kg . 

The predictions for the slope of the symmetry energy Lq 
and the nucleon Dirac mass given in table Q] for the 
two models are consistent with experimental values lying 
in the range 60-90 MeV [16] for the slope and in the range 
0.7-0.9A/at for the Dirac masses. 

Hypcrons are included as in Ref. [10] : the coupling con- 
stants arc written as gf n — w{ n g m where is the weight 
factor for the meson m coupled to the baryon /. For 
the mesons cj, p, tt and 5, the meson-hyperon couplings 
are deduced from the meson-nucleon couplings based on 
SU(6) symmetry imposing w£ = (1 + s^/3), where 
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Table 2: Weight factor w{ for the hyperons A, S and S. 



Model 


A 


S 




MC1-H 


0.61 


0.18 


0.33 


MC1-HF 


1.00 


0.50 


0.55 


MC2-HF 


1.00 


0.51 


0.52 


MC3-HF 


0.91 


0.38 


0.48 


MC4-HF 


0.89 


0.39 


0.47 


QMC-H 


0.62 


0.35 


0.32 


QMC-HF1 


0.71 


0.51 


0.37 


QMC-HF2 


0.64 


0.41 


0.33 


QMC-HF3 


0.70 


0.33 


0.37 



is the strangeness of the baryon /. For the isovector p 
or pseudo-vector n mesons, this symmetry is accounted 
for by the isospin and Gamow-Teller operators. Only the 
coupling constants of the scalar s meson are allowed to 
vary around the values imposed by SU (6) symmetry. In 
QMC700, the weights w{ and w{ are adjusted in the bag 
model and scale with the bag radius |10j . In the chiral 
model and in the extensions of QMC in table [IJ they are 
adjusted to the non-relativistic potentials of hyperons at 
saturation, 

V NM ) = Zi + j^Zi + %^, (4) 

where is the single-particle energy and Y> s and Sg are 
the scalar and time-component self-energies Vat,r(/) is 
fixed to be -30 MeV for / = A [T7], 30 MeV for / = £ [TO] 
and -18 MeV for / = H 0]. While the value for the A is 
quite certain, that for the other hyperons are still under 
debate. The weights w{ are given in table [2] while we 
choose w{ = (1 + s-^/3) for all baryons /. 

Being partly related to the underlying quark structure 
of baryons the chiral and QMC models are good candi- 
dates for exploration of the properties of dense core of neu- 
tron stars. In table [3] we explore different models for the 
EoS with and without hyperons and their predictions for 
the maximum mass M max and associated radius R(M max ). 
The first three EoS are based only on nucleons. For the 
following three models in table the EoS is based on 
non-interacting hyperons. The maximum mass for these 
models is strongly reduced compared to models without 
hyperons. Interacting hyperons arc taken into account in 
the next 10 models in table [31 For the five chiral models 
the maximum mass is increased compared to the models 
with free hyperons, but it remains smaller than that pre- 
dicted by the models without hyperons. There is a small 
variation of about 0.05 M© between the different HF mod- 
els depending on the extension of the chiral model under 
consideration. 

In the last four lines of table [3] we present variants of 



Table 3: Different models corresponding to different approxi- 
mations for nucleons (second column) and for hyperons (third 
column). The last three columns correspond to the maximum 
mass M max , its associated radius -R max and the corresponding 
baryonic density n max - 



EoS 


Nucleon 


Hyperon 


Ail ax 


fl(M max ) 


max 




N 


Y 


(Af©) 


(km) 


(in no) 


MC1-H/N 


MC1-H 


No 


1.65 


10.6 


7.8 


MC1-HF/N 


MC1-HF 


No 


1.87 


10.8 


7.2 


MC3-HF/N 


MC3-HF 


No 


1.85 


10.7 


7.3 


MC1-H/NY FG 


MC1-H 


Fermi Gas 


1.02 


11.9 


5.6 


MCI-HF/NY^g 


MC1-HF 


Fermi Gas 


1.12 


12.0 


5. 1 


MC3-HF/NY FG 


MC3-HF 


Fermi Gas 


1.13 


11.9 


5.5 


MC1-H/NY 


MC1-H 


MC1-H 


1.45 


11.9 


6.5 


MC1-HF/NY 


MC1-HF 


MC1-HF 


1.58 


12.2 


5.6 


MC2-HF/NY 


MC2-HF 


MC2-HF 


1.60 


12.2 


5.5 


MC3-HF/NY 


MC3-HF 


MC3-HF 


1.55 


12.0 


5.6 


MC4-HF/NY 


MC4-HF 


MC4-HF 


1.55 


12.0 


5.6 


QMC700/NY 


QMC700 


QMC700 


2.00 


12.5 


5.1 


QMC-H/NY 


QMC-H 


QMC-H 


1.68 


12.2 


5.5 


QMC-HF1/NY 


QMC-HF1 


QMC-HF1 


1.97 


12.3 


5.2 


QMC-HF2/NY 


QMC-HF2 


QMC-HF2 


1.97 


13.3 


1.6 


QMC-HF3/NY 


QMC-HF3 


QMC-HF3 


1.66 


12.3 


5.3 



the QMC model whose parameters are given in table [T] 
The original QMC700/NY model predicts a maximum 
mass of 2.00 M with hyperons [10]. The model QMC- 
H/NY is a consistent Hartrcc model where the s and uj 
coupling constants are readjusted to reproduce saturation 
properties (no, Bq and V/v.r(/))- The maximum mass in 
the model QMC-H/NY is lower than that of the origi- 
nal QMC700/NY since the incompressibility is lower. In 
the model QMC-HF1/NY the Hartree and Fock terms are 
consistently treated in the mean-field potential and the 
same readjustment procedure as in QMC-H/NY is fol- 
lowed. The Fock contribution to the total energy is how- 
ever treated non-relativistically as in the QMC700/NY 
model. Including the complete Fock contribution in the 
energy, we obtain the fully relativistic QMC-HF2/NY 
model. The maximum mass is not altered by the fully 
relativistic treatment of the Fock term. Finally, we re- 
duced the incompressibility by increasing the parameter 
C from 0.51 to 1 (see table [T| and the maximum mass 
decreases from 1.97 to 1.66 M©. From this analysis, we 
conclude that the large mass obtained in the QMC700/NY 
model |10j is mainly due to an overestimated and incom- 
pressibility. 

We show in figure [T] the mass-radius relation for the 
various EoS presented in this Letter. Dotted lines show 
results for the models without hyperons while the solid 
lines stand for the models including interacting hyperons. 
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Fig. 1: 
EoS. 



Radius 



Mass of neutron stars versus the radius for selected 



As expected the EoS with hyperons reach lower maxi- 
mum masses than the models without hyperons, except 
the model QMC700 due to its large and unrealistic incom- 
pressibility. The maximum mass obtained with hypernu- 
clear EoS does not change much for different models. This 
is due to a self-regulating compensation effect between the 
softening of the EoS and the onset of the hyperonic degree 
of freedom [7|5] . The maximum mass of about 1.6 M© pre- 
dicted by the relativistic approach is slightly larger than 
that of 1.4 Mq found in Ref. [8]. The difference of 0.2 Mq 
is mainly related to the larger incompressibility of the rel- 
ativistic models. It is also observed that the radii for the 
Hartree models are smaller than that of the Hartree-Fock 
ones, and that the predicted radii of the relativistic models 
and that of Ref. [8] are comparable. 

The abundances for nucleons and hyperons are shown 
for various EoS in figure [2j with free hyperons (MC1- 
H/NY FG ) and with interacting hyperons (MC1-H/NY 
and MC1-HF/NY). The order of appearance of hyperons 
is different for these models and depends strongly on the 
corresponding interactions. In the model MCl-H/NYj? 
with free hyperons, the S~ is non interacting and appears 
first. On the other hand, in the models MC1-H/NY and 
MC1-HF/NY, the E~ interaction is repulsive while the 5~ 
interaction is attractive, therefore the latter tends to pre- 
vent the £~ from appearing. In the model MC1-H/NY, 
the A appears first followed closely in density by the S~, 
then at a larger density by the 3°, while for the model 
MC1-HF/NY the 5~ and 5° appear first and push the 
other hyperons to larger densities. 

We have presented in this Letter different EoS for hy- 
pcron matter derived from both chiral and QMC models. 
These models are thcrmodynamically consistent and treat 
the Hartree and Fock terms on equal footing. We have 
found that the different chiral parameterizations give a 
maximum neutron star mass that does not exceed 1.6 Mq 
for HF models with K w 270 MeV. Correcting the 
QMC700 model [TU] by including the complete relativistic 
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Fig. 2: Particules abundances for the various models of table[3] 
MC1-H/NY FG (top), MC1-H/NY (middle) and MC1-HF/NY 
(bottom). 



contribution to the Fock term, a correct incompressibility 
and the experimental constraints for the hyperon mean- 
fields, we found a decrease of the predicted maximum mass 
from 2.00 down to 1.66 Mq. We additionally checked that 
the effect of rotation with the period of 3.15 ms (as for PSR 
J1614-2230) does not increase the maximum mass by more 
than 0.02 Mq. We therefore conclude that it is difficult 
to reconcile any of the model presented in this work with 
the observed mass of PSR J1614-2230 and the empirical 
knowledge of saturation properties of nuclear matter. Un- 
less an unexpected property of the hyperon interaction is 
missing in the present model, as well as in the BHF models 
of Refs. [6"r [8l2"rj] . our results tends to exclude hypernuclear 
matter to be present in the core of massive neutron stars. 
It reinforces alternative models such as for instance the 
model of deconfined quark matter, for which experimen- 
tal constraints are almost inexistent. We can therefore 
conclude that the equation of state of dense matter in the 
core of neutron star is still not even qualitatively under- 
stood. 
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